Dust in the wind: Crystalline silicates, corundum and periclase in 

PG 2112+059 



F. Markwick-Kemper/'^ S. C. Gallagher,^ D. C. Hines*^ & J. Bouwman^ 



ABSTRACT 

We have determined the mineralogical composition of dust in the Broad Ab- 
sorption Line (BAL) quasar PG 2112+059 using mid-infrared spectroscopy ob- 
tained with the Spitzer Space Telescope. From spectral fitting of the sohd state 
features, we find evidence for Mg-rich amorphous sihcates with ohvine stoichiom- 
etry, as well as the first detection of corundum (AI2O3) and periclase (MgO) in 
quasars. This mixed composition provides the first direct evidence for a clumpy 
density structure of the grain forming region. The silicates in total encompass 
(56.5 ±1.4)% of the identified dust mass, while corundum takes up (38 ±3) wt.%. 
Depending on the choice of continuum, a range of mass fractions is observed for 
periclase ranging from (2.7 ± 1.7)% in the most conservative case to (9 ± 2)% in 
a less constrained continuum. In addition, we identify a feature at 11.2 fim as 
the crystalline silicate forsterite, with only a minor contribution from polycyclic 
aromatic hydrocarbons. The (5 ± 3)% crystalline silicate fraction requires high 
temperatures such as those found in the immediate quasar environment in order 
to counteract rapid destruction from cosmic rays. 

Subject headings: galaxies: active — quasars: emission lines — quasars: individ- 
ual (PG 2112+059) — ISM: dust, extinction 
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Introduction 



In the local Unive rse, Asympto tic Giant Branch (AGB) stars are believed to dominate 



dust production (e.g., iDraind l2003l ). AGB stars are a late phase of stellar evolution for 
stars with initial masses of M=1-8Mq, which develop outflows where the low temperatures 
and high densities become optimal for dust condensation. However, a significant amount of 
dust originating from AGB st ars does not build up with in ~ 1 Gyr of the birth of the first 
generation of low-mass stars (IMorgan fc Edmunds! l2003l ). Nevertheless, large quantities of 
dust are clearly present at these early times: quasar host galaxies at 2; ~ 6 show evidence 
for 10^~^ Mq of dust he ated by star format i on in their host g alaxies as seen in submm 
and far-infrared emission JPriddev et allboosl : iBeelen et al.l 120061 ). Furthermore, extinction 
curves f or a 2: = 6.2 quasar and GRB 050904 {z = 6.29) are not consistent with those lo- 



cally (IMaiolino et al.ll2004l: IStratta et al 



20071 ). Another mechanism for dust production is 



therefore required, and lElvis et al.l (|2002| ) explored the possibility that quasar winds might 
provide environments suitable for efficient dust formation. They determined that the tem- 
peratures and pressures in these regions could reach the values found in cool, dus t-producing 



stars. The quasar winds are predicted to produce dust masses up to 10'' Mq ( lElvis et al. 



2OO2I ). implying that the quasar wind may in some cases only account for part of the dust 
producti on, and supernoyae ha ve been suggested as an alternative source of dust in high-z 
galaxies (ISugerman et al.l 120061 ). 



Quasar outflows are most obviously manifested in Broad Absorption Line (BAL) quasars 



This population, approxirn ately 20% of optically selected type 1 quasar samples (e.g.. lHewett &: Foltz 
2OO3I : iReichard et al.ll2003l ). is notable for broad, blueshifted absorption evident in common 
ultraviolet resonance transitions such as C IV, Lya, and O VI. These P Cygni-type features 
arise because the observer is looking through an outflowing wind. BAL quasars are thus a 
natural population to consider when investigating the grain properties of dust in the quasar 
environment. In this letter, we present a detailed analysis of the mid-infrared spectrum 
of the luminous BAL quasar, PG 2112+059, in order to determine its dust composition. 
PG 2112+059, an IRAS source, is mid-infrared bright, and one of the most luminous low- 
redshift {z = 0.466) Bright Quasar Survey objects with My = —26.9 (IIo= 70kms~^ Mpc~^, 
flu = 0.3, and = 0.7 are assumed throughout). HST spectra reveale d broad, shallow 
C IV absorption, and it has been well-studied in the UV and X-ray (e.g., 
2004h . 



Gallagher et al. 



2. Observations and Data Reduction 



PG 2112+059 was observed on May 25, 2005, using the low resolution modules of 
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the InfraRed Spectrograph (IRS; iHouck et al.ll2004l ) on board Spitzer ( IWerner et al.l 12004 ). 
under AOR key 10949376. We performed 1 cycle of 240 sec. integrations with the short low 
modules (5.2-14.5 /im) and 2 cycles of 120 sec. with the long low modules (14.0-38.0 /xm). 

We app lied the data reduction , background sub t ractio n and flux ca l ibrati on methods 
described in lBouwman et al.l (120061 ): iBouwman et al.l (120071 ): ISwain et al.l (120071 ) to pipeline 
version S14 of the data. The resulting absolute flux is accurate to ~5%, while the relative 
fluxes are accurate to ~1.7%. The resulting combined spectrum is shown in Fig. [H 



3. Analysis and Results 



The spectrum of PG 2112+059 exhibits the b road ernission bands at ~ 10 and ~ 18 
/im generally ascribed to amorphous silicates (e.g. iDraind 120031 ). In addition, the [Ne ll] 
transition at 12.8 fim is visible, as well as substructure in the solid state emission at ~ 11.2 
/im (Fig. [T]). 



3.1. Continuum Determination 



as 



Complementary to the IRS spectrum, we added mid-IR photometry from the literature 
compiled by NEE0 (see Figure [1]). In wavelength regions in com mon, we found goo d 



agreement between (all observed-frame w avelengths) 25/im IRAS d ata ([Sanders et al 



12 and 7.5/im ISO data JHaas et al.lboooh . and 10/im Palomar data (INeugebauer et al.l 



1989), 



1987h . 



To determine the continuum underlying the dust emission features, we fit a power-law model 
{Fi, = Fy^ov") using the photometric and spectroscopic data (and uncertainties) from 5-8/im 
and 24.5-40/im(rest-frame). As shown in Figured! this provided a good fit to the continuum 
with a = —0.617 ± 0.004. We note that there is no evidence from quasar cor aposite spectral 



energy distributions for an inflection under the dust emission features (see [Richards et al. 



20061 ). and a power-law model fits the luminous quasar composite in this spectral regime 
quite well. Because the continuum is not well constrained at wavelengths longer than the 
broad 18/xm emission feature, we performed a second fit including the 13-14/im IRS data at 
the dip between the 10 and 18/im emission features. This pulled the continuum up to the 
value of the 60/im /i?y4 5" photometric datapoint plus uncertainty for a = —0.674 ± 0.003. 



^NASA/IPAC Extragalactic Database; 



http : / / nedwww . ipac . caltech . edu 
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3.2. Spectral Feature Analysis 



Assuming the emission region is optically thin in the infrared, the continuum-divided 
spectrum (Fig. [2]) yields the resulting opacity of the dust. As at each wavelength, features 
and continuu m contributions are in equ al amounts arising from the same temperature com- 
ponents (e.g. ISiebenmorgen et al.ll2005l ). the temperature dependence has been eliminated 
by dividing the the spectra by the continua. The c ontribution of sy nchrotron emission to 
the continuum at these wavelengths can be ignored (IHaas et al.lll998l ). 



We fitted the resulting opacity curves with laboratory spectroscopy of minerals com- 
monly observed in AGB stars, assuming that the resulting opacity is a linear combination 
of the individual mineral opacities. We also subtracted a spline-fit continuum between our 
fit boundaries - 8 and 25 fim - of each mineral spectrum, to provide the same baseline. 
For each continuum divided spectrum, we ran 512 models with varying dust compositions 
(amorphous and crystalline olivine, AI2O3, MgO, polycyclic a romatic hydr ocarbons (PAHs); 
based on prior observations and model calculations, see e.g. iDraind l2003l ) and grain shape 
(spherical or non-spherical; iBohren fc Huffman! 1 198 31 ). to find the relative mass fractions in 
each component. The fit was performed using minimization of each continuum-divided 
spectrum, excluding the 12.8 /im [Nell] line. 

Fig. [2] shows the best fits for both continuum-divided spectra, with Ax^ values of 2.8 for 
continuum 1 and 2.6 for continuum 2. If we allow a tolerance of a factor 2 for Ax^, we find 
classes of dust compositions containing 44 and 124 good fits; 30 and 47 of which require the 
presence of PAHs over the 8-25 fim range, for the first and second continuum respectively. 

The essential dus t components in t his class of fits ar e amorphous olivine - MgFeSi04 



( Dorschner et al. 1995) and/orMg^SiO^ 



1999I ). corundum (Al90q: lBegemann et al. 



ager et al.ll2003l) crystallin e forsterite (IKoike et al. 
I997I) and periclase fMgO: lHofmeister et al.ll2003h . 



For the 20 /xm excess in the spectrum not explained by the silicates, we find that MgO in 
non-spherical grains (represented by a common distribution of ellipsoids) provides a good 
match. Both spherical and non-spherical grains give good results for the amorphous olivines 
and corundum. Including grains with radii a ^ 1 /xm greatly increased the Ax^ of the fits, 
and are therefore not likely to contribute to the observed emission features. The fit results 
are summarized in Table [H 

The dust mixture is very similar for both continuum-divided spectra, with the notable 
exception of the MgO mass fraction. In case of the second continuum, where the 13-14 
micron range was included in the continuum fit, the MgO contribution is suppressed, but 
when the 13-14 micron range is -more realistically- treated as dust feature emission, and is 
not used to constrain the continuum, we find that ~ 9 wt.% of the dust may be in the form 
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of MgO. 



To measure the contri bution of PAHs to the emission, we used an average interstellar 
profile for the PAHs from iHony et al.l (120011 ) . Although they contribute to the 11.2 fim 
feature, the PAHs do not reproduce its shape well, and additional opacity from forsterite is 
needed. We find that PAHs are responsible for ~30% and ~15% for the first and second 
continuum respectively. This is consistent with the appearance of the PG 2112+059 spectrum 
in the 5-9 fim range (Fig. [2]). The 11.2 /xm resonance in PAHs is likely due to neutral 
mole cules, while the features in th e 5-9 fim range are predominantly carried by ionized PAHs 



(e.g. Ivan Diedenhoven et al.ll2004l ). A range of a factor ^^6 in /11.2//6.2 around the interstellar 
value is observed in Galactic environments covering more than 4 orders of magnitude in Go 
( IHony et al.ll200ll ). In case Gq is significantly stronger than the prevalent field in the Galaxy, 
I11.2/ Iq.2 decreases, further limiting the amount of PAHs contributing to the 11.2 /xm feature. 



4. Dust composition 



The presence of silicates in local AGN has been known for rn any years through res- 
onant absorption at 9.7 and 18 /xm (e.g., Ilmanishi fc Uend I2OOOI ). while more recently. 



silicate features have also been detected in emission towards several quas ars and AGN 
dSiebenmorgen et al.l[2005l : Isturm et allbood : iHao et al.l[2005l : Ishi et al.lEooel l 



While the spectral appearance of PG 2112+059 in the infrared clearly requires silicates, 
it is interesting that oxides appear to be present as well. As a source of opacity, corundum 
and periclase cause the 18 /xm feature to become stronger relative to the 9.7 /im resonance, 
while at t he same tirne, the y do fill in the opacity trough between the two silicate resonances 
(see also IJager et al.l l2003l ). Corundum also causes a red wing on the 9.7 /xm resonance, 
which fits well with the feature as observed in PG 2112+059. This is the first reported 
detection of periclase and corundum as dust components in a quasar. 



4.1. Corundum 



The presenc e of Al20,s is riot sur prising; it is commonly found in the stellar winds of 



AGB stars (e.g. ISloan &: Pried Il998l ). In thermodynamic equilibrium, corundum (Al^O^ 
is ex pected to condense at ~ 1500 K, well above the stability limit for silicates (ITielens 
19901 ). In a cooling gas newly condensed AI2O3 grains will gradually be covered by silicate 
layers as the temperature drops, provided the density of the gas does not decrease too 
rapidly. In cases where the gas density has become too low for significant dust growth once 
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the gas temperature has reached the stabihty hmits for sihcates, a freeze-out of the dust 
condensation sequence occurs, and AI2O3 remains the dominant dust component. Indeed, 
AGB star s with low mass loss rates (10~^ - 10~^ M (7^ yr~^) contain significant fractions 
of AI2O3 ( Blommaert et al. 2006 : Heras fc Honyl 2005 ). and a decr ease in AI2O3 conten t is 
observed with increasing mass loss rate for AGB stars in the LMC ( iDijkstra et al.ll2005l ). 

Both t he luminosity an d prese nce of silicates in emission in PG 2112+059 (rather than 
absorption; ISmith et al.l e.g., 120071 ) a rgue for geomet rically associating these grains with the 
quasar rather than the host galaxy (jShi et al.l 120061 ). The presence of both corundum and 
olivines in the spectrum indicate that there is a highly inhomogeneous - clumpy - density 
str ucture in the dust fo rming regions as expected fro m the broad-band spectral modeling 
of iNenkova et al.l (120021 ) and the theoretical work of lElitzur fc Shlosmanl (120061 ). As the 
gas travels outwards, certain pockets apparently maintain high enough density throughout 
the cooling process that silicates can be formed on top of the corundum, whereas in lower 
density regions of the wind the dust condensation sequence gets truncated after the formation 

of AI2O3. 



4.2. Periclase 



Periclase (MgO) has its stability limits ~50 K below the stability limits of Mg-rich 
olivines (~1100 K), but just above the evaporation temperature o f Fe - the temperature 



below which Fe^^ can be incorporated in the amorphous silicates (IGail fc Sedlmayd Il999 



Ferrarotti fc Gailll2003l ). In the thermodynamic equilibrium of a slowly cooling gas, Mg2Si04 
would form first, co nsuming all available ga s-phase Mg in the process, thus prohibiting the 
formation of MgO JPerrarotti fc Gaillboosh . The presence of solid MgO in PG 2112+059 
indicates that the dust forming gas cooled rapidly to a temperature below these stability 
limits, allowing for the formation of MgO and Fe-containing silicates at the same time as 
the formation of Mg2Si04. 



4.3. Amorphous silicates 

The amorphous silicates in PG 2112+059 have the stoichiometry of a Mg-rich olivine. 
By using both Mg2Si04 and MgFeSi04 to fit the spectrum, we have determined that the 
Fe-content of the amorphous silicates is low: Mgi.95Feo.o5Si04, although using only Mg2Si04 
or MgFeSi04 gives good results, too. The silicate grains are found to be sub-micron-sized 
(at 10 /xm the opacity of silicates is the same for all grain sizes < 0.5 /xm), but the data did 
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not allow us to distinguish between spherical or non-spherical grains. These results compare 
well to the sil icates in the Galact ic interstellar medium,which are composed of spherical 
olivine grains (jKemper et al.ll2004l ). or even more general, of non-spherical Mg-rich olivines 
dMin et al.ll2007l l 

Attempts to determine the silicate properties in active galaxies include a fit to the 9.7 
/im absorption feature seen in NGC 1068, which shows evidence for the presence of Al- 
bearing species, either in the f orm of aluminurn -silicate s or corundu r n, giv ing rise to a red 
wing to the silicate resonance (jjaffe et al.l 120041 ) , while ISturm et al.l (120051 ) notice that the 
9.7 and 18 /xm emission features seen in NGC 3998 have a different spectral appearance 
than the silicates in the Galactic ISM. In particular, the 18 /xm feature is very strong with 
respect to the 9.7 /im feature, and the peak positions o f the resonances deviate as well. We 
suggest that in NGC 3998, as well as in some quasars (jshi et al. 2006 ). the broad 'silicate' 
features might be more accurately modeled with a mixture of silicates, corundum and MgO, 
not unlike PC 2112+059. 



4.4. Crystalline silicates 



A fraction of (5±3)% by weight of the silicates is found to be crystalline. Evidence of 
crystalline silicates is not reported in any other silicate emission spectra seen in quasars or 
AGN, though a feature seen at 11.2 /im in some objects is identifie d with PAHs, supporte d 



by the detection of additional PAH resonances at 6.2 and 7.7 (jSchweitzer et al. 



20061). 



We only ascribe a small fraction of the emission around 11.2 fim to PAHs (see Sect. 13. 2p . 
Cryst alline silicates are no t detected in the interstellar medium of our own Galaxy (< 2.2 
wt.%; iKemper et al.ll2005l ). due to radiation damage caused by cosmic rays on a time scale 
of 40 Myr. In quasars, the times scales may be even shorter if the cosmic ray flux is higher, 
and crystalline silicates need to be reformed continuously to explain the observed abundance. 
The formation of crystalline silicates requires higher densities than amorphous silicates. A 
quasar wind origin may explain the presence of cryst alline silicates in a relatively small 
fraction (12 out 77) of ultraluminous infrared galaxies (ISpoon et al.ll2006l ). 



5. Conclusion 

For the first time, the composition of the dust in a quasar has been determined, albeit 
limited to the components with resonances in the infrared. While the origin of the dust could 
lie in stellar ejecta, the dust properties are consistent with their formation in the quasar 
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wind itself. The dust in PG 2112+059 clearly bears similarity to dust in other astrophysical 
environments, i.e. in the properties of the amorphous silicates, but the presence of large 
amounts of MgO and AI2O3 sets it apart from the composition of interstellar dust in the 
local universe. 

Assuming the dust is formed in the quasar wind, the co-existence of the highly refractory 
corundum, the crystalline silicates and and the less refractory MgO and amorphous silicates 
indicates that the wind of PG 2112+059 has an inhomogeneous temperature and density 
structure. 

This work is based on observations made with the Spitzer Space Telescope, which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by NASA through an award issued 
by JPL/Caltech. 

Spitzer (IRS) 
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Table 1. Best A^^ fits to the two continuum-divided spectra of PG 2112+059 (Fig. E]). 
Tlie first column shows the dust fraction considered, and the following columns show the 
mass fractions for both continuum fits and the average. The standard deviation on the 
results is determined from the spread of the mass fractions observed in those fits within the 
tolerance of a factor 2. The bottom line shows the average composition of the 

amorphous silicate dust. 





continuum 1 


continuum 2 


average 




wt.% 


wt.% 


wt.% 


silicates/total 


(54 ± 1) 


(59 ± 1) 


(56.5 ±1.4) 


corundum/total 


(37 ±2) 


(39 ±2) 


(38 ± 3) 


MgO/total 


(9 ±2) 


(2.7 ± 1.7) 


(5.9 ±2.6) 


crystalline / silicates 


(4 ±2) 


(6 ±2) 


(5 ±3) 


Mg2Si04/am. silicates 


(70 ± 44) 


(49 ± 40) 




MgFeSi04/am. silicates 


(30 ± 44) 


(51 ±40) 




am. sil. composition 


Mgi.9Feo.iSi04 


Mgi,8Feo.2Si04 


Mgi.85Feo.i5Si04 
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Fig. 1. — {Top panel) Mid-infrared spectrum of PG2112-I-059; the abscissas are rest-frame 
values. The diamonds are photometric data from the fol lowing sources (from left to right, 
given in observed -frame units): I RAS 60 and 25 um ({Sanders et aL 1989), ISO 12 um 



feaas et al.lboooh . Palo mar 10 /im (iNeugebauer et al 
Palomar 3.5/xm (two observations; INeugebauer et al. 



■9871) ISO 7.5 /xm flHaas et al.ll2000h . 
19871 ). Shown is the IRS spectrum 



(with error bars), and the dashed and thin solid curves are two power-law continuum models 
{Bottom panel) The IRS spectrum minus the dashed continuum fit (continuum 1; thick solid 
histogram) and thin solid continuum fit ( continuum 2; thin histogram) shown in the top 
panel. The crosses indicate the frequencies of the data used in the dashed continuum fit; 
for the second continuum fit, the IRS data from 13-14 fim ([2.30-2.14] x 10^^ Hz) were also 
included. 
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Fig. 2. — Best 8 — 25 /xm fits, from our suite of good fits, for the continuum-divided spectrum 
of PG 2112+059 (offset by -1) using tlie first and second clioice of continuum (top and 
bottom panel, respectively). The thick solid line is the total fit to the data (histogram, with 
error bars) composed of amorphous olivine (dotted line), forsterite (dashed line), corundum 
(AI2O3, dash-dotted line), MgO (dash-triple-dotted line), and the mean interstellar PAH 
spectrum (solid line; only in panel 1). The inset in panel 1 shows the strength of PAH 
features expected in the 5.5-8.5 /im range for a radiation field comparable to that of the 
Galaxy. 



